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Abstract. Graphene is a gapless semiconductor in which conduction and valence 
band wavefunctions differ only in the phase difference between their projections onto 
the two sublattices of the material's two-dimensional honeycomb crystal structure. We 
explain why this circumstance creates openings for broken symmetry states, including 
antiferromagnetic states in monolayer and bilayer graphene and exciton condensates 
in double-layer graphene, that are momentum space analogs of the real-space order 
common in systems with strong local interactions. We discuss some similarities among, 
and some differences between, these three broken symmetry states. 
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1. Introduction 

Graphene[l,j 2j is a two-dimensional crystal consisting entirely of carbon atoms. Its 
honeycomb lattice is stabilized primarily by strong planar sp^ bonds, leaving one more 
weakly bonded vr-orbital per carbon atom available for metallic conduction. Because of 
the two atoms per unit cell in a honeycomb lattice the tt electrons form two bands, one 
of which is occupied in a neutral sheet. Because the vr-bands cross at the inequivalent 
honeycomb lattice Brillouin-zone corners, i.e. at the K and BZ corner points, the 
gap between the occupied tt valence band and the unoccupied tt conduction band 
vanishes. Most graphene properties of interest to condensed-matter physicists, for 
example transport and optical properties, directly involve only the tt electron orbit als 
that are close to the band crossing points. These are accurately described over an energy 
interval several eV in width by a • p Hamiltonian that has the form of a massless Dirac 
equation: 

In Eq.Q Vq is the velocity of band electrons at the band-crossing (Dirac) point which 
can be related to band Hamiltonian vr-electron hopping amplitudes, and a are Pauli 
matrices which act on sublattice labels s{s^). Eq.Q applies near the valley K Dirac 
point; the corresponding equation for valley is obtained by letting kx —kx- When 
the sublattice degree-of-freedom is viewed as a pseudospin, we see from Eq.Q that the 
band eigenstates are chiral. This is in the conduction band of valley K pseudospin 
is parallel to the momentum, while in the valence band pseudospin is antiparallel to 
momentum. 

Because only one of the two band states is occupied at each momentum, the 
many-body ground state can be continuously deformed relative to the non-interacting 
ground state without breaking translational symmetry simply by rotating the pseudospin 
direction at each momentum. When interactions are neglected the valence band is full, 
implying that the pseudospin direction is always opposite to the direction of momentum. 
Since the interaction energy is minimized when all pseudospins are parallel, as we 
will discuss explicitly below, there is tension between band energy minimization and 
interaction energy minimization in monolayer graphene. A similar tension arises in two 
other graphene based two-dimensional electron systems - graphene bilayers|^ which 
consist of two Bernal stacked graphene layers and graphene double-layers [4j which 
consist of two layers separated by an insulating tunnel barrier. In this article, we 
will discuss the competition between interaction and band energies in all three systems 
using a common mean-field language which enables us to highlight similarities and point 
out differences. We will discuss the potential broken symmetry states using a mean- 
field-theory in which the ground state is determined by minimizing the total energy of 
single-Slater-determinant many-body states varying the pseudospin direction at each 
momentum k. By performing a stability analysis for the resulting energy functional we 
conclude that broken symmetries states can occur in monolayer graphene if interactions 
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are sufficiently strong, and that they occur in bilayer and double-layer graphene for 
interactions of any strength. Our mean-field theory treats the gapped [5l |6l [7] and 
nematic states|^i9j that have been discussed for bilayer graphene on an equal footing. 

A parallel can be drawn between the vr-orbital states in these three graphene systems 
and the electronic states of Mott insulators. In Mott insulators strong interactions 
project most of the many-electron wavefunction onto a subspace with one atom per unit 
cell. The only degree-of-freedom that is available in this subspace is the spin state at each 
lattice position. Unless interactions between spin orientations on different lattice sites 
are frustrated, the ground state is normally close to the classical ground state in which 
spin-orientations are fixed on each site and chosen to minimize the expectation value of 
the effective spin Hamiltonian. Quantum fiuctuations of spin orientations play only a 
quantitative role. For the graphene states discussed here the pseudospin energy function 
that is fixed by energy minimization plays a role similar to the spin distribution in an 
insulator, but is a function of momentum rather than position. The graphene system 
band Hamiltonians, which differ essentially in the three cases, act like momentum- 
dependent pseudospin fields. Because the band Hamiltonians reduce symmetries, energy 
minimization does not in all cases imply broken symmetry ground states. Because of 
the absence of a gap in the graphene case, important quantum fiuctuations occur in 
both pseudospin orientations and in the occupation numbers of momentum states, but 
we expect that their role is also only quantitative and we do not discuss them at length 
in this paper. The broken symmetry states which can occur in these graphene systems 
are unusual from several points of view. For a given spin and valley, the ordered states 
of both monolayer and bilayer graphene have large quasiparticle Berry curvature and 
spontaneous quantized anomalous Hall effects [lOl [HI Ej. The potential broken 

symmetry states of double-layer graphene have spontaneous interlayer phase coherence, 
which leads to the suite of phenomena connected with counter-fiow superfiuidity [151 ES] 
when the layers are separately contacted. 

Our paper is organized as follows. In Section 2 we formulate the mean-field 
theory in the three graphene-based two-dimensional electron systems in terms of energy 
minimization with respect to fc-dependent pseudospin orientation. In Sections 3, 4, and 5 
we discuss a stability analysis of the energy functional at the band state configuration for 
monolayer, bilayer, and double-layer systems respectively and comment on the character 
of the broken-symmetry states that result. Finally in Section 6 we point to important 
similarities and differences between the three cases. 

2. Pseudospins in Monolayer, Bilayer, and Double-Layer Graphene 

We consider a variational single-Slater-determinant wavefunction with a single 
pseudospin state occupied at each momentum fc, and minimize the expectation value of 
the Hamiltonian with respect to the pseudospin orientations h^. In this approximation 
the four spin-valley fiavors interact only through their contribution to the electrostatic 
Hartree energy. Since we will consider only states that are locally electrically neutral. 
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we will ignore the Hartree energy for the moment and return to it later when we discuss 
the role of spin- valley flavors. 




I^> + \B) \B) 

Figure 1. Representation of a pseudospin state by a unit vector r{6^(j)). Pseudospins 
directed along the ±z direction represent the states labeled \A) and \B). For monolayer 
graphene these states are confined to alternate sublattices, while for bilayer and double 
layer graphene they are confined to alternate layers. General pseudospin states are 
coherent linear combinations of \A) and \B)^ with the azimuthal pseudospin angle 
(/) specifying the phase difference between \A) and \B) amplitudes and the polar 
pseudospin angle specifying the \A)-\B) polarization. A state with = 7t/2 has equal 
\A) and \B) weights, whereas the states with ^ = and = tt equal \A) and \B) 
respectively. 



It will be convenient to express the band Hamiltonian as an effective magnetic fleld 
that acts on the pseudospin degree-of-freedom by writing 

'Hband= Yl ^\s'^h'^s's)%s- (2) 

k,s' ,s 

In this language the three cases we discuss are distinguished by their pseudospin effective 
magnetic flelds: 

h^^^^ = hvok [cos(0^) X + sin(0^) y] , 

hf^^ = -{h^e/2m') [cos(2(/),^) X + sm{2(t>^) y] , 

hf^^ = hvQ{k-kF)z, (3) 

for monolayer, bilayer, and double-layer cases respectively. Here 0^ is the angular 
orientation of the two-dimensional k momentum and k is its magnitude. Note that the 
band Hamiltonians are off-diagonal in pseudospin index in the monolayer and bilayer 
cases, and diagonal in the double-layer case. The pseudospin labels in Eq.([2]) refer to 
the sublattice index in the monolayer graphene case and to layer index in the bilayer 
and double-layer cases. The form we have chosen for the bilayer Hamiltonian applies 
only at energies smaller than the interlayer hopping energy 71 and is due [3] to virtual 
hopping between the two low-energy bilayer vr-orbital sites, which are located in different 
layers, via two higher-energy vr-orbital sites that are not explicitly retained. Because the 
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pseudospin labels refer to position in all three cases, the electron-electron interaction 
Hamiltonian is diagonal in pseudospin at each vertex in all three cases: 

^int = ^ E E 4+.>4-.>'Cp>'Ca:,. [v^iQ) + V-{q)a:,a:,,] (4) 

k,p,q 

where V^{q) = {Vs ± Vb)/2 and Vs^b are the momentum-space interactions between 
electrons on the same and different sublattices (or layers) respectively. In the monolayer 
case the interactions are pseudospin independent (Vs = V^)^ whereas in the bilayer and 
double-layer cases they are pseudospin dependent because interlayer interactions are 
weaker than intralayer interactions. 

The dependence of the band energy on the momentum-dependent pseudospin 
configuration is easy to evaluate. From the form of the pseudospinor for a state aligned 
in a particular pseudospin direction (see Fig.[l]), 

^-{:$*)' 

it is easy to show that 

{h\a\h) = n , (6) 

and that 

\n){n\ = . (7) 

In Eq.([5]), 9 and (j) are the polar and azimuthal angles for the pseudospin orientation, 

i.e. h = (sin 9 cos 0, sin 9 sin 0, cos 9). For the case of pseudospin-dependent interactions 

the following identity is useful: 

. 1 — a^n^ — a^ny + a^n^ , , 

a |n)(n|a = ^ . (8) 

Taking the expectation values of the band and interaction Hamiltonians we find that 

k 

and that 

^int[%] = ^ E [(1 + Vs{k -p} + {nln^^ + ny.) V^k - p}] . (10) 

k,p 

The interaction energy is an exchange contribution which sets the momentum transfer q 
in Eq.li]) to p—k. As explicitly shown in Eq.( 10), the exchange energy contribution from 
any pair of momentum is lowered when their pseudospins are made more parallel. The 
band energy, on the other hand, is minimized when the pseudospin direction is opposite 
to the direction of at every k and hence strongly pseudospin dependent. In the 
following sections we will discuss the pseudospin orientation function which minimizes 
the total energy E'tot = ^'band + ^int for monolayer, bilayer, and double layer cases. 
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3. Spin-Density- Wave States in Monolayer Graphene 

We find it useful to expand the total energy functional to leading order around its band 
theory value ^i^g = This consideration follows similar lines in the three cases of 

interest. To leading order we can preserve normalization by writing 

n^k = -k{^-lfy')+^k (11) 

where is a two-dimensional vector perpendicular to h^. For monolayer graphene 
the band pseudospin orientation n^g = —k = — (^cos 0^, sin 0^, 0^ and 5 has z and 

azimuthal (0) components along z = (0,0,1) and = z x k = (^sin 0^, — cos ^g, 0^ 
directions respectively We find that 

E = Eo + ^EESIKlpl (12) 
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Figure 2. Pseudospin orientations and quasiparticle spectra near the Dirac point for 
gapless and gapped monolayer graphene. The upper row illustrates conduction band 
pseudospin orientations for gapless (a-left) and gapped (b-right) cases; the valence 
band pseudospins are opposite in direction at each momentum. The arrows indicate 
the magnitude and direction of the x — y plane pseudospin projection while the color 
indicates the z direction pseudospin projection which is zero in the gapless unbroken 
symmetry state. In monolayer graphene pseudospins rotate around the z axis by 
27r when a point in momentum space encloses the Dirac point /c = 0. In the broken 
symmetry state the self-energy contribution to the pseudospin field has a component in 
the z direction which does not vanishes for A: ^ 0, opening up a gap in the quasiparticle 
spectrum. 
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where 



where 



= + - Mk -p)4>r 4'f -Vs(k- P) s<,A (13) 

(14) 

The quantity in the above equations, the change in interaction energy associated 
with switching a single band state from valence to conduction band pseudospin 
orientations, can be identified as the exchange contribution to the self-energy of the band 
state. adds to the energy difference between conduction and valence band states 
because of the energy cost of reversing pseudospin orientation at a single-momenta, 
keeping all other pseudospins fixed. For the monolayer graphene case we can let 
Vb ^Vs 80 that 

^k = ^Yl ^(P- ^) ' (15) 

p 

which vanishes for a 5-function interaction model because of the angular average over 
the direction of but grows with \k\ for the realistic Coulomb interaction case: 
d 

hvok ln(A/fc) (16) 

where a = e^/ehv^ is graphene's fine structure constant, A ^ 1/a is the Dirac 
model ultraviolet cutoff where a is graphene's lattice constant, and e is the graphene 
sheet's effective dielectric constant which depends on the substrate used to support the 
sheet. Full Brillouin zone Hartree-Fock theory calculations suggest that the most 
appropriate value for A is ^ 30/a. The self-energy term captures the physics of a 
theoretically anticiapted [18l [I9[ EOj EH |22] logarithmic interaction enhancement of the 
energy difference between conduction and valence band quasiparticles, that has now 
been confirmed experimentally [23] . This effect that is normally described in terms of 
an interaction enhanced quasiparticle velocity at momenta near the Dirac points. 

Depending on the dielectric environment of a monolayer graphene sheet, the fine 
structure constant value can vary between a ^ 0.5 and a ^ 2[24j. It is sometimes 
claimed that perturbative and mean-field treatment of electron-electron interaction 
effects, like the one discussed here, cannot be trusted because the coupling constant 
is not small. A more reliable way to judge the adequacy of these approximations 
is to compare with experiments. In the case of monolayer graphene application of 
this type of criteria usually argues for the opposite conclusion, namely that mean-field 
theory is reliable, provided that the electron-electron interaction is properly screened 
when carriers are present. The approximation in which electronic self-energies are 
approximated at leading order in dynamically screened Coulomb interactions, variously 
known as the random phase approximation (RPA) or the GW approximation, agrees 
very well with, in particular, photoemission experiments j25j in both neutral and charged 
monolayer graphene. For neutral graphene the inverse quasiparticle lifetime [261 |27j in 



Order in Graphene 



8 



neutral graphene is ^ 10% of the quasiparticle energy at typical values of a, and this 
ratio may provide a better characterization of interaction strength. 

Although some details of graphene interaction physics are sensitive to the long 
range of Coulomb interactions [18l EHl [29l [30], most importantly perhaps the velocity 
enhancement mentioned above, we are able to make a number of valuable points by 
considering a short-range interaction model in which the momentum dependence of Vs 
and 1/d are neglected. Note that because we are neglecting inter- valley scattering we are 
still imagining that the interaction range is large compared to the atomic length scale, 
and that there is therefore no direct relationship between this approximation and using 
a lattice Hubbard model 03l El SSI SSj . Setting Vs = Vb ^ U we find that S vanishes, 
that 



i^l'f ^ 2AhvokS^,^. - U^^, ■ 4>f, (17) 



and that 



Kg.^2AhvokSj:.-U. (18) 

In-plane and out-of-plane pseudospin reorientation instabilities are indicated by 
vanishing eigenvalues for Kpt and K^^^ respectively. We search for a zero eigenvalue by 
solving 

\EkI;;s^^^o. (19) 

p 

These homogeneous equations are solved by setting Si C co^{(f)^ — x)/^'^ofc and 
5^ C/hvok where C is an arbitrary constant and x an arbitrary angle. We obtain the 
conditions 

2co.(4-x)^t//(g, '°-'^'--y^---^' (20) 

for /C^'* and 

for K^'^. Converting the integral over momentum p into an integral over the energy of 
quasiparticles with energy hvop^ it follows that instability occurs in K"^'"^ when [/z/* > 4, 
whereas for K^'^ it occurs when [/z/* > 2. In these equations z/* = W/{27r?i^VQ) is the 
Dirac-model density-of-states at the model's ultraviolet energy cutoflF scale W ^ hvo/a; 
the integrand of the energy integral is constant because of a cancellation between 
the hvop factor in the denominator and the quasiparticle density-of-states, which is 
proportional to energy. Interactions are less effective in reducing the energy cost of 
in-plane distortions than out-of-plane z distortions because of the angle-dependence 
of the band state pseudospins. 

There are two important points to make about these stability criteria results. First 
of all, we see that if interactions are strong, pseudospins are more likely to tilt toward 
the ±z directions, rather than to alter their orientations in the x — y plane. The state 
produced by this pseudospin distortion has a higher electron density on one sublattice 
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than the other, hence it is a density-wave state when looked at from a microscopic point 
of view. In ignoring the Hartree energy, we have imphcitly assumed that the signs of 
the density- waves are opposite for opposite vaUeys or for opposite spins, with the latter 
possibility being more likely as discuss later. The expected state is therefore a spin- 
density- wave state rather than a charge- density- wave state. If the instability involved 
distortion of the in-plane pseudospin direction, rather than tilting out of the plane, it 
would yield a state with spontaneous anisotropy, characterized by the angle in which 
the magnitude of the quasiparticle self-energy depends on momentum direction. The 
mean-field-theory instability analysis therefore suggests that density-wave instabilities 
occur before anisotropy instabilities. Secondly, we note that the instability criterion 
involves z/*, the density-of-states at the band width energy scale. We can conclude from 
this observation that the presence or absence of a broken symmetry state depends on 
atomic length scale physics beyond that captured by the k • p Dirac model. [131 113 HH] 
When solved with an unscreened Coulomb interaction [5], mean-field theory predicts 
that the instability occurs in monolayer graphene at a ^ 1.3; screening and other higher- 
order corrections will shift [3T] the instability - likely to larger values of a beyond those 
that can be reached in monolayer graphene even when it is suspended so that interactions 
are not reduced by dielectric screening. Degrees of freedom not included in the vr-band 
only model of graphene, and portions of the Brillouin-zone far from the the corner Dirac 
points are also likely to play an important role. The instability can also be ruled out 
experimentally with nearly complete confidence because pseudospin orientations in the z 
direction open gaps in the quasiparticle spectrum. It can be established experimentally 
the gaps, if present, cannot be larger than ~ 10~^ eV compared to the natural energy 
scale of graphene which is ^ 10 eV. [32j 



4. Antiferromagnetic States in Bilayer Graphene 

Several consequential distinctions can be drawn between the monolayer and bilayer 
cases. First of all the two-band k • p for bilayer model provides a good description 
only at energies that are small compared to the interlayer tunneling energy 71 ^ 0.4 
eV, whereas the corresponding monolayer two-band model applies up to energies 
^ 2 eV, much closer to the full band width. Secondly, both the magnitude and 
direction of pseudospin effective magnetic fields have diflFerent behaviors since ^^^^^ 
varies quadratically rather than linearly with k and has an orientation angle that is 
twice the momentum orientation angle. The stability analysis for bilayers parallels 
that for monolayers precisely once these diflFerences are recognized. For bilayers the 
band and in-plane-distortion pseudospin directions are n.t = (cos(2(/)g),sin(2(/)g),0) 

and 0g = n^g X z = (^sin(2(/)g), — cos(2(/)g), 0^. The more rapid variation of pseudospin 
direction with momentum direction eliminates the logarithmic divergence of the velocity 
enhancement at small k found in the monolayer Coulomb interaction case, but still [22j 
leaves a substantial interaction-induced velocity enhancement. For the short-range 
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interaction model, the stability matrices for cj) and z distortions are 



(22) 



and 



27.2 



where 0^ • now equals co8{2(/)^ — 20^). Setting (7cos(2(/)p' — x)/(^^P^/2m*) 

and 5|> ^ C / {Ti?p^ /2rrf) we can solve for the interaction strength at which the 
smallest eigenvalue approaches zero. As in the monolayer case we find that the ±z 
distortions, which in this case correspond to moving charge between layers, occur 
at weaker interaction strengths. The instability criteria [6j in the bilayer case are 
Uu^^\n{W/EF) > 4 for K^'^ and Uiy^^ln{W/ Ep) > 2 for K'^'. In this case 
^BL ^ m*/(27r/i^) is the energy-independent band electron density-of-states. The 
1/E quasiparticle energy factor is not canceled by an increasing density-of-states, as 
it was in the monolayer case, and the interaction contribution to the stability eigenvalue 




Figure 3. Pseudospin orientations and quasiparticle dispersions near the Dirac 
point in bilayer graphene. The upper row illustrates the conduction band pseudospin 
orientations; the valence band pseudospin orientation is opposite at each momentum. 
In bilayer graphene pseudospins rotate by an angle of An around the z direction when 
a point in momentum space encloses the ^ = Dirac point. The figures (a) and 
(b) represent respectively gapless and gapped broken symmetry states. As in the 
monolayer case all pseudospins are in the xy plane in the ungapped state. In the 
gapful case the ±i out of plane pseudospin components are non-zero and indicate 
sublattice polarization. In the low energy model for bilayers the sublattice index labels 
layer. 
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integral has a logarithmic infrared divergence which we have cut off by assuming that 
conduction band states up to energy have been occupied, Pauli-blocking pseudospin 
polarization. For the conclusion is that the density-wave instability occurs 

before the anisotropy distortion, and that it will occur for arbitrarily weak interactions 

mm- 

5. Exciton Condensates in Double-Layer Graphene 




Figure 4. Pseudospin orientations and quasiparticle dispersions in double layer 
graphene near the Dirac point. In the two-band pseudospin model the two bands 
furthest from the Fermi energy are not accounted for explicitly. In double layer 
graphene tunneling between layers is negligible so that the pseudospin label is 
equivalent to a layer label. The band pseudospin direction changes abruptly between 
up and down directions along the common Fermi surface where the conduction band of 
the high density layer and the valence band of the low density layer cross. In presence 
of interactions a gap is opened and pseudospin directions rotate gradually between up 
and down directions. The figures (a) and (b) represent respectively the single particle 
bands crossing at a Fermi circle, and the gapped phase with interlayer coherence. A 
gap opens in presence of arbitrarily weak interactions as in the case of bilayer graphene. 
Any in plane pseudospin component introduces interlayer coherence and reduces the 
total energy of the system. In our illustration we have chosen the interlayer coherence 
pseudospin component to point in the x direction. 

Now we turn to the case in which two graphene layers, one containing electrons in 
the conduction band and the other containing an equal density of holes in the valence 
band, are coupled by repulsive Coulomb interactions. If we ignore the completely full and 
completely empty energetically remote bands, we can view the double-layer graphene 
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system using the same pseudospin language that we used for the monolayer and bilayer 
graphene cases. In the absence of interactions the conduction band states of the n- 
layer, which we associate with pseudospin up, are occupied inside a Fermi circle and 
the valence band states of the p-layer, which we associate with pseudospin down, are 
occupied outside of this circle (see Fig. 4). As in the monolayer and bilayer cases, one 
pseudospin state is occupied at each momentum. Because the band pseudospins are 
oriented in the ±z directions, rather than in the x — y plane, the pseudospin distortions 
that can potentially lower the energy are in the x and y direction distortions, rather 
than d) and z distortions. We find that 



= = [^Mhvoik - k^) + T.^)5^,^. -Vo{k- p)\ , (24) 



where 



nf r 



= ^ E Vsip- k) (1 + nl^). (25) 

P 

(We take the band energy to be the quasiparticle energy in the absence of carriers in 
either layer and include a self-energy contribution from the full valence band of the 
n-type layer.) The sudden change in pseudospin orientation at the Fermi circle has 
a cost in exchange energy which can be mitigated by rotating pseudospins into the 
X — y plane, which corresponds to establishing coherence between layers spontaneously. 
Because K^'^ = X^'^, the pseudospin rotation can occur with the same gain in energy at 
any azimuthal angle. The fact that the energy is independent of the interlayer phase (z.e, 
the pseudospin azimuthal angle) implies that this broken symmetry state supports super 
currents that fiow in opposite directions in opposite layers [U [IGj. As in the monolayer 
and bilayer cases, instability is indicated by a vanishing K^'^ eigenvalue. The instability 
condition can be found by solving 

^1:3-0 (26) 

p 

with 5| C/hvo\k — For a 5-function interaction the self-energy term simply 

shifts the relationship between Fermi energy and density and plays no role. The 
instability criterion is therefore Uiy^^ln{2EF/6) > 2 where we have chosen 2Ef as 
an ultraviolet cutoff, S is an infrared cut-off, and = Ep / {2tiTi?v'^) is the constant 
density-of-states of the double-layer Dirac model for energies between and Ep. As 
in the bilayer case, an instability occurs for arbitrarily weak interactions. Although 
this conclusion is universally accepted by researchers who have examined this possible 
ordered state [JJ [33l |34l [35l [36] estimates of the size of the consequent energy gap vary 
widely because of the difficulty of accounting accurately for the infiuence of carrier 
screening. 



6. Discussion 



Graphene two-dimensional electron systems (2DES) are remarkable for several different 
reasons. The fact that they are truly two-dimensional on an atomic length scale 
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elevates 2DES physics from the low-temperature world to the room-temperature world. 
Furthermore, they are accurately described by very simple models over very wide energy 
ranges and yet have electronic properties that can be qualitatively altered simply by 
stacking [3, [371 EHl ESI SQj them in different arrangements, and by adjusting external 
gate voltages. In this article we have discussed the properties of three different graphene 
2DES's using a simple mean- field-theory pseudospin language and specializing to the 
case of electrically neutral systems. 

The basic building block of all graphene 2DES's is the isolated monolayer, which 
is described by a massless Dirac k • p Hamiltonian over a wide energy range. The 
Dirac model has chiral quasiparticles, and in the graphene case the chirality refers to 
the relationship between k • p momentum and the direction of a pseudospin associated 
with the sublattice degree-of-freedom of graphene's honeycomb lattice [5l E]. In 
neutral graphene each momentum is singly occupied on average. In our mean- field- 
theory approach we do not account for quantum fluctuations in these momentum 
occupation numbers, but allow energy to be minimized with respect to the pseudospin 
orientation at each momentum. Using this approach we find that in monolayer graphene 
strong interactions can lead to a broken symmetry state in which the pseudospin 
rotates from the x — y plane toward the ±z directions, breaking inversion symmetry 
and opening up a gap in the quasiparticle excitation spectrum. This conclusion has 
been reached previously |42l |43l [44[ |45l [46l [IE] by several researchers, sometimes using 
more sophisticated theoretical approaches which attempt quantitative estimates of the 
required interaction strengths. In k - p mean-field theory the broken symmetry state 
consists of a density-wave state with more charge on one sublattice than the other 
within each spin-valley flavor, but no overall charge density variation. Since exchange 
interactions beyond the k level favor[13j states with the same sublattice polarization 
on both valleys, the strong interaction state is likely a spin-density-wave state. As 
it happens, it appears to be clear from experiment that this broken symmetry state 
does not occur^23l [25j in monolayer graphene, even when suspended. This property 
is generally consistent with theoretical expectations and is perhaps unfortunate from 
the point of view of researchers interested in many-body phenomena. Happily the very 
closely related bilayer and double-layer graphene systems that were not anticipated 
prior to the experimental emergence of the graphene field, are more likely to have 
broken symmetry states and may save the day for interesting many-body phenomena in 
graphene based 2DES's. 

Bilayer graphene is described, over a more limited energy range however, by a 
similar pseudospin k • p model with quadratic rather than linear dispersion; in bilayer 
graphene the quasiparticle velocity vanishes with momentum as in a conventional 2DES, 
but the pseudospin direction still depends on the direction of momentum. In pseudospin 
language, it is clear that the bilayer is more susceptible to instabilities to broken 
symmetry states because the cost in band energy of rotating the pseudospins toward 
the ±z direction is smaller at small momentum. Indeed we find that instabilities occur 
for arbitrarily weak interactions O [6]. In the bilayer case the spin- density- wave state 
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has opposite pseudospin orientations for opposite spins [T2l [131 HI]- An alternate state 
in which the pseudospin orientation is rotated within the x — y plane to increase the 
degree of alignment leads to an anisotropic state, and gains less exchange energy for 
a given interaction energy cost. There is indeed a great deal of evidence from recent 
experiments [4^ [49] that a broken symmetry state does occur in bilayer graphene, but 
the character of the state is not yet completely settled since some experimenters find 
evidence for a gapless anisotropic state [47j and others ^Ql [5l] evidence for a gapped 
isotropic stated [6l [71 [I3]. On the theoretical side different researchers have also reached 
different conclusions concerning the character of these states, with some researchers [8^ 9j 
concluding that the broken symmetry state should be anisotropic. 

In mean-field theory the gapped broken symmetry state has lower energy than the 
anisotropic broken symmetry state because all band pseudospinors can be tilted to the 
±z direction to reduce the probability of finding electrons with parallel pseudospins. The 
efficacy of in-plane pseudospin distortions is reduced because their pseudospin-directions 
(j)^ depend on the momentum orientation angle 0^. To us, the conclusion that the 
weak-coupling broken symmetry state will be gapped appears to be unavoidable unless 
somehow overturned, in a way which has not yet been clearly articulated, by inter-fiavor 
correlations. In drawing conclusions from the renormalization group calculationsj6l[8l[9] 
which attempt to go beyond the mean-field theory considerations described here, it is 
important to realize that, because short-range interactions within a valley act only 
between opposite pseudospins, the pseudospin dependence of the corresponding fiowing 
interaction has no significance in a many-fermion Hilbert space. 

Bilayer graphene differs from monolayer graphene mainly because of its weaker 
dependence of band energy on pseudospin direction at small momentum. This difference 
is sufficient to lead to states with broken pseudospin symmetry. Double-layer graphene 
differs in a more qualitative way because not only the interaction energy, but also the 
band energy, is diagonal in the z component of pseudospin, i.e. in layer index. The band 
Hamiltonian in this case has a sudden change in the sign of the z direction pseudospin 
orientation. This momentum space domain wall has a large interaction energy cost 
which can be mitigated by rotating the pseudospins near the Fermi surface out of the 
z direction with a common azimuthal angle. In this way, the momentum dependence 
of the pseudospin rotates smoothly between z and —z directions, and the interaction 
energy is lowered. Like bilayer graphene, double-layer graphene has a broken symmetry 
for arbitrarily weak interactions. In both cases the size of the gap is difficult to estimate 
quantitatively [U [33l [3ll [35l [36] and likely to be overestimated by mean-field-theory The 
quasiparticle density-of-states remain finite at low energies in the double-layer graphene 
case because the pseudospin field vanishes along a line in momentum space while in the 
bilayer case because it varies quadratically with wavevector near = 0, but the end 
result is essentially the same. 

In two-dimensional electron systems, the quantized Hall conductance can be 
expressed [52] in terms of an integral of Berry curvature over momentum space. Using 
ideas from topology it is possible to show that the Hall conductivity is always e'^/h times 
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an integer valued topological invariant known as Chern number. For two-band models, 
momentum states can always be described using a pseudospin-1/2 language like the one 
used in this paper. For this case the quantized Hall conductivity is especially simple 
to visualize geometrically since it is equal to the number of times the unit sphere of 
pseudospin directions is covered upon integrating over momentum space. The broken 
symmetry states of bilayer graphene carry a positive or negative unit of Hall conductivity 
because they cover either the north-pole or the south pole twice when the pseudospin 
at A; = points to the north or south pole. (In the double-layer graphene case, on 
the other hand, the pseudospin is confined to the interface between a plane and the 
pseudospin unit sphere, i.e. to a line on the unit sphere that encompasses zero area.) 
The inversion symmetry breaking states in bilayer graphene can therefore be viewed as 
spontaneous quantum Hall states [HI [121 ESj. In the spin-density- wave state the Hall 
contributions from different valleys cancel, but states with a non-zero Hall conductivity 
can be stabilized by going to a non-zero total carrier density in weak magnetic fields. 

In this article we have discussed the possibility of interaction driven broken 
symmetries in three different graphene based two-dimensional electron systems, 
monolayer graphene systems and two layer graphene systems that are stacked in two 
different ways. The three cases we have discussed are most easily addressed theoretically, 
but not by any means the ones that are most likely to have strong interaction effects 
and broken symmetries. ABC stacked multilavers [37l [381 139] . for example, tend to have 
even smaller separations between conduction and valence bands at small but are 
complicated by competing electronic structure details. Recent advances in techniques 
for preparing samples in which disorder plays an inessential role appear to be bringing 
us close to clear experimental conclusions as to the strength and character of broken 
symmetries in bilayers and trilayers. These developments will, no doubt, reveal some 
surprises that present some focused challenges to theory. 
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